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1.  General  Program  Description 

The  program  outlined  in  our  original  proposal  is  directed  to  the  development  of 
constructive  methods  for  certain  constrained  optimization  problems  arising  in  the 
design  and  control  of  electromagnetic  fields  and  in  the  identification  of  scattering 
objects.  The  problems  to  be  addressed  fall  into  three  categories:  (i)  the  design  of 
antennas  with  optimal  radiation  characteristics  measured  in  terms  of  directivity; 
(ii)  the  control  of  the  electromagnetic  scattering  characteristics  of  an  object,  in 
particular  the  minimization  of  its  radar  cross  section,  by  the  choice  of  material 
properties;  and  ( iii )  the  determination  of  the  shape  of  scattering  objects  with  various 
electromagnetic  properties  from  scattered  field  data.  As  originally  envisciged.  the 
main  thrust  of  the  program  is  toward  the  development  of  constructive  methods 
based  on  the  use  of  complete  families  of  solutions  of  the  time-harmonic  Maxwell 
equations  in  the  infinite  domain  exterior  to  the  radiating  or  scattering  body.  During 
the  course  of  the  work  an  increasing  amount  of  attention  has  been  devoted  to  the 
use  of  iterative  methods  for  the  solution  of  various  direct  and  inverse  problems.  The 
continued  investigation  and  development  of  these  methods  and  their  application  in 
parameter  identification  has  Isecome  a  significant  part  of  the  program. 

2.  Research  Activity 

Considerable  attention  during  the  course  of  the  program  has  been  diri-cted 
toward  the  application  of  optimization  techniques  to  the  problem  of  determining 
currents  on  conformal  anteniicis  which  optimize  various  functionals  of  the  radiated 
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field  such  as  directivity  and  signal-to-noise  ratio.  Physical  limitations  on  the  design 
of  radiating  structures  impose  certain  restrictions.  One  consequence  of  which  is  that 
the  optimization  problems  are  subject  to  constraints. 

Our  work  (see  e.g.  [4],  [26j)  has  demonstrated  the  fetisibility  of  employing 
complete  families  of  solutions  in  these  constreiined  optimization  problems  in  which 
the  boundary  data  serve  as  control  parameters.  Indeed,  not  only  have  we  done  a 
basic  analysis  of  this  class  of  problems,  but  also  have  developed  computer  codes  for 
particular  problems. 

We  were  invited  t )  prepare  a  review  of  these  methods  for  a  special  issue  of  the 
Proceedings  of  the  lEE.'H  and  a  manuscript  will  appear  shortly  which  demonstrates, 
both  theoretically  and  computationally,  how  a  variety  of  antenna  performance  cri¬ 
teria  may  serve  as  cost  functionals  for  the  general  approach.  .\  copy  of  our  review 
paper  is  attached.  Dissertation  work  in  this  area  has  been  carried  out  by  two  stu¬ 
dents.  S.  Fast  and  T.  Bazow.  The  Ph.D.  thesis  of  Fast  An  Opinnization  Mnihod 
for  Solving  a  Rodiation  Dirextion  Problem,  included  codes  for  the  three  ilimensioiial 
electromagnetic  problem  of  directing  a  signal  into  a  prescribed  three  dimensional 
sector  by  appropriate  choice  of  surface  current.  T.  Bazow.  who  is  currently  finishing 
his  work,  is  studying  the  more  complicated  problem  of  usmg  surface  impedance  as 
the  control  parameter  for  the  problem  of  controlling  the  scattered  field  due  to  an 
incoming  plane  wave.  Unlike  the  acoustic  and  electromagnetic  cases  studied  earlier, 
the  problem  is  not  reducible  to  a  generalized  eigenvalue  problem.  The  cejnstructive 
procedure  for  finding  optimal  impedances  is  therefore  considerably  more  difficult 
and  is  again  based  on  the  use  of  complete  families  of  solutions. 

Such  optimal  radiation  problems,  if  meaningfully  formulated,  fall  naturally  into 
the  claj5S  of  constrained  optimization  problems.  These  constraints  typically  arise  in 


two  ways:  inherent  constraints  due  to  physical  limitations  of  the  system  e.g.  bounds 
on  input  power,  and  artificial  constraints  whose  introduction  is  motivated  by  the 
desire  to  control  unwanted  effects.  For  example,  it  has  long  been  recognized  that 
the  narrow  focusing  of  the  main  beam  of  an  antenna  has  the  concomitant  (and 
undesirable)  effect  of  increasing  near  field  power.  Indeed  one  typically  wishes  to 
both  focus  the  wave  beam  and  minimize  the  power  stored  in  the  wave  field.  Thus  a 
typical  problem  that  arises  in  antenna  design  is  that  of  dealing  with  several  possibly 
conflicting  goals. 

Such  problems  may  be  cast  in  the  framework  of  multi-criteria  optimization. 
.Apparently,  such  a  formulation  has  not  been  made  for  problems  in  the  design  of 
radiating  structures  although  the  techniques  are  known  and  have  been  applied  to 
a  variety  of  engineering  design  problems  especially  in  mechanical  engineering.  We 
have  prepared  a  manuscript  [41]  proposing  the  use  of  these  techniques  and  giving 
examples  of  numerical  calculations  for  a  multi-criteria  problem  related  to  that  of 
maximizing  the  signal-to-noise  ratio. 


We  have  also  been  concerned  with  establishing  existence,  uniqueness  and  con¬ 
tinuous  dependence  results  for  certain  "non-classical”  bovmdary  value  prolDlems 
in  both  the  acoustic  and  electromagnetic  problems  of  time-harmonic  scattering. 
Thus,  in  the  acoustic  case  we  have  considered  the  so-called  resistive  and  conductive 
"boundary”  conditions  to  model  scattering  from  thin  shell-like  structures.  This 
work  has  been  extended  to  electromagnetic  problems  where  our  formulation  in  fact 
generalized  the  usual  conditions  for  the  transmission  problem.  In  addition  to  well- 
posedness  we  have  also  considered  properties  of  the  set  of  far  field  patterns.  This 
work  is  reported  in  [7],  [13],  [14],  [29],  [34]. 

Each  of  these  topics  of  boundary,  or  more  accurately  generalized  transmission. 
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conditions  is  intended  to  more  accurately  model  the  behavior  of  scattering  surfaces 
which  are  penetrable  to  electromagnetic  radiation  and  serve  as  better  approxima¬ 
tions  to  the  physical  situation  of  scattering  from  control  bodies.  The  use  of  such 
coatings,  modelled  by  boundary  data  such  as  resistivity,  for  the  control  of  scattered 
radiation  represents  an  area  of  mathematical  research  comparatively  new  and  of 
potentially  significant  application.  The  work  done  under  this  grant  is  preparatory 
to  the  study  of  optimizing  scattering  characteristics,  e.g.  the  minimization  of  radar 
cross  section,  by  choice  of  impedance,  conductivity  or  resistivity. 

We  have  also  devoted  considerable  effort  to  inverse  problems,  both  of  shape 
reconstruction  and  profile  inversion.  With  respect  to  the  former  problem,  we  have 
cleveloped  a  constructive  method  using  complete  families  of  radiating  solutions. 
These  feunilies  are  used  to  formulate  the  inverse  problem  of  choosing  the  shape 
whose  scattering  characteristics  best  approximate  the  measured  far  field  data  as  a 
penalized  optimization  problem  whose  penalization  term  involves  only  the  defect  in 
satisfying  the  boundary  conditions.  In  collaboration  with  B.  Kok.  we  have  developed 
a  working  computer  code  for  the  problem  of  reconstructing  the  shape  of  a  perfectly 
conducting  cylinder  when  illuminated  by  an  E-polarized  incident  electromagnetic 
wave.  Work  is  being  completed  by  X.  .Jiang  on  a  three-dimen.sional  code  for  the 
TM  case  and  a  manuscript  describing  the  algorithm  and  the  computational  results 
is  being  prepared.  Results  of  this  aspect  of  the  work  have  been  reported  in  [4].  [6]. 
1201,  (21|,  [25|,  126], 

The  importance  of  a  good  initial  guess  in  the  shape  reconstruction  algorithm 
has  become  apparent  in  the  course  of  the  numerical  work,  particularly  in  the  full 
three  dimensional  code.  One  way  to  obtain  a  rea.sonable  starting  surface  is  to 
carry  out  the  inversion  in  a  simpler  setting.  This  can  be  done  at  low  frequencies 


4 


and  work  wat;  done  for  reconstruction  problems  in  the  low  frequency  regime.  ;xs 
well  as  for  solutions  to  the  direct  problem.  Not  only  has  an  easily  programmable 
algorithm  been  developed  to  find  the  best  ellipsoidal  approximation  to  an  unknown 
scattering  shape,  but.  remarkably,  it  was  found  that  Kelvin  inversion  could  be  used 
for  solutions  of  the  Helmholtz  equation,  contrary  to  accepted  belief,  thus  enlarging 
the  class  of  surfaces  for  which  solutions  could  be  found.  These  methods  can  then 
be  used  to  refine  the  estimate  of  the  scattering  surface  used  as  a  starting  point  m 
the  inverse  problem.  This  work  on  Kelvin  inversion  has  appeared  in  two  papers  by 
G.  Dassios  and  R.  E.  Kleinman  [5],  [9]. 

The  algorithm  for  reconstruction  of  the  best  approximating  ellipsoid,  originally 
formulated  for  the  case  where  measurements  of  the  entire  far  field  are  available  [1], 
has  been  recently  modified  for  the  case  in  which  only  backscattered  data  is  available, 
and  a  computer  code  has  been  written  and  tested  by  X.  Jiang,  using  a  forward 
code  based  on  the  exact  low  frequency  result  for  ellipsoids.  We  intend  to  continue 
testing  with  experimental  data  supplied  by  colleagues  at  the  University  of  Michigan 
Radiation  Laboratory  who  are  collaborating  under  an  NSF  grant.  Work  is  now 
in  progress  to  compile  this  algorithm  with  the  original  inverse  scattering  code  to 
provide  a  good  initial  surface  for  the  reconstruction  algorithm. 

Work  continued  on  the  second  class  of  inverse  problems  involving  parameter 
estimation.  In  particular  we  considered  the  problem  of  reconstructing  an  unknown 
complex  index  of  refia<.,tion  from  known  scattered  field  data.  This  problem  of  elec¬ 
tromagnetic  tomography  or  profile  inversion  has  received  considerable  attention  in 
recent  years.  In  its  usual  formulation  it  is  a  constrained  nonlinear  optimization  prolj- 
lem  in  which  both  the  field  and  the  index  of  refraction  are  unknown  in  the  scattering 
object.  Moreover  the  field  depends  on  the  index  of  refraction  in  a  highly  nonlinear 
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way.  Most  attacks  on  this  problem  embody  two  principles:  first  a  linearization  of 
the  nonlineax  dependence  and  secondly  a  regularization  of  the  ill-posed  optimiza¬ 
tion  problem.  The  process  is  an  iterative  one  which  requires  the  direct  solution  of 
a  scattering  problem  at  each  step  in  order  to  determine  the  updated  field. 

Our  approach  was  motivated  by  our  work  on  iterative  solutions  of  the  integral 
equations  based  on  error  minimization.  This  work  resulted  in  a  uniform  framework 
which  encompassed  a  number  of  methods  -  successive  over  relaxation,  conjugate 
gradient  and  Krylov  methods.  Applications  of  these  various  methods  to  the  integral 
equations  of  scattering  was  also  considered.  This  work  is  described  in  [2],  [3],  [8], 
[10],  [11],  [12],  [23],  [24],  [27],  [30],  [31],  [33],  [35],  [39]. 

One  of  the  simplest  methods  in  this  class  is  successive  over  relaxation  where 
at  each  step  the  relaxation  parameter  is  determined  so  as  to  minimize  residual  er¬ 
ror.  We  incorporated  this  idea  in  the  inverse  problem  by  adding  the  constraint, 
the  state  equation,  as  a  regularizer.  We  then  proposed  a  modified  gradient  method 
to  solve  the  optimization  problem.  Both  the  field  and  the  profile  were  updated  at 
each  step  by  adding  a  relaxation  parameter  times  an  updating  function  (direction). 
The  direction  for  the  field  was  taken  to  be  the  residual  error  in  the  state  equation, 
exactly  as  in  the  over-rela.xation  solution  of  the  direct  problem,  while  the  direction 
for  the  index  of  refraction  was  taken  to  be  the  gradient  of  the  error  in  matching 
the  measured  data.  The  relaxation  parameters  were  found  simultaneously  by  mini¬ 
mizing  the  regularized  functional.  This  required  the  solution  of  nonlinear  algebraic 
equations  for  the  unknown  parameters  arising  from  the  necessary  conditions  for  a 
minimum.  The  novel  features  of  this  approach  are  the  elimination  of  the  need  to 
solve  a  direct  problem  at  each  step  while  retaining  the  nonlinear  nature  of  the  cost 
functional.  This  work  was  carried  out  in  collaboration  with  P.  M.  van  den  Berg  of 
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Delft  University  of  Technology  in  the  Netherlands.  Preliminary  numerical  results 
aje  very  promising  and  have  been  reported  in  [16]  and  [37].  Recent  results  on  two 
dimensional  objects  indicate  that  the  method  is  capable  of  reconstructing  discon¬ 
tinuous  profiles  even  with  noisy  data.  The  following  figures  show  cm  original  profile 
and  the  reconstructions  after  64  iterations  where  the  data  consisted  of  the  scattered 
field  at  30  points  equally  spaced  on  a  circle  of  radius  3A  produced  by  a  line  source 
located  successively  at  each  of  the  30  locations. 


original  reconstruction  reconstruction 

no  noise  10%  noi.se 


manuscript  describing  these  results  is  under  preparation. 
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